Structural insights into the activation mechanisms of human HtrA serine proteases.

Human HtrA1-4 proteins belong to the HtrA family of evolutionarily conserved serine proteases and function as important modulators of many physiological processes, including maintenance of mitochondrial homeostasis, cell signaling and apoptosis. Disturbances in their action are linked to severe diseases, including oncogenesis and neurodegeneration. The HtrA1-4 proteins share structural and functional features of other members of the HtrA protein family, however there are several significant differences in structural architecture and mechanisms of action which makes each of them unique. Our goal is to present recent studies regarding human HtrAs. We focus on their physiological functions, structure and regulation, and describe current models of activation mechanisms. Knowledge of molecular basis of the human HtrAs' action is a subject of great interest; it is crucial for understanding their relevance in cellular physiology and pathogenesis as well as for using them as targets in future therapies of diseases such as neurodegenerative disorders and cancer.

Distinct 3D Architecture and Dynamics of the Human HtrA2(Omi) Protease and Its Mutated Variants.

HtrA2(Omi) protease controls protein quality in mitochondria and plays a major role in apoptosis. Its HtrA2S306A mutant (with the catalytic serine routinely disabled for an X-ray study to avoid self-degradation) is a homotrimer whose subunits contain the serine protease domain (PD) and the regulatory PDZ domain. In the inactive state, a tight interdomain interface limits penetration of both PDZ-activating ligands and PD substrates into their respective target sites. We successfully crystalized HtrA2V226K/S306A, whose active counterpart HtrA2V226K has had higher proteolytic activity, suggesting higher propensity to opening the PD-PDZ interface than that of the wild type HtrA2. Yet, the crystal structure revealed the HtrA2V226K/S306A architecture typical of the inactive protein. To get a consistent interpretation of crystallographic data in the light of kinetic results, we employed molecular dynamics (MD). V325D inactivating mutant was used as a reference. Our simulations demonstrated that upon binding of a specific peptide ligand NH2-GWTMFWV-COOH, the PDZ domains open more dynamically in the wild type protease compared to the V226K mutant, whereas the movement is not observed in the V325D mutant. The movement relies on a PDZ vs. PD rotation which opens the PD-PDZ interface in a lid-like (budding flower-like in trimer) fashion. The noncovalent hinges A and B are provided by two clusters of interfacing residues, harboring V325D and V226K in the C- and N-terminal PD barrels, respectively. The opening of the subunit interfaces progresses in a sequential manner during the 50 ns MD simulation. In the systems without the ligand only minor PDZ shifts relative to PD are observed, but the interface does not open. Further activation-associated events, e.g. PDZ-L3 positional swap seen in any active HtrA protein (vs. HtrA2), were not observed. In summary, this study provides hints on the mechanism of activation of wtHtrA2, the dynamics of the inactive HtrA2V325D, but does not allow to explain an increased activity of HtrA2V226K.

Potent antidiuretic agonists, deamino-vasopressin and desmopressin, and their inverso analogs: NMR structure and interactions with micellar and liposomic models of cell membrane.

Deamination of vasopressin (AVP) enhances its antidiuretic activity. Moreover, introduction of D-Arg8 instead of its L enantiomer in deamino-vasopressin (dAVP) results in an extremely potent and selective antidiuretic agonist - desmopressin (dDAVP). In this study we describe the synthesis, pharmacological properties and structures of these two potent antidiuretic agonists, and their inverso analogs. The structures of the peptides are studied in micellar and liposomic models of cell membrane using CD spectroscopy. Additionally, three-dimensional structures in mixed anionic-zwitterionic micelles are obtained using NMR spectroscopy supported by molecular dynamics simulations. Our conformational studies have shown that desmopressin in a membrane mimicking environment adopts one of the characteristic for vasopressin-like peptides ß-turn - in position 3,4. Furthermore, dDAVP shows the tendency to create a ß-turn in the Cys6-Gly9 C-tail, considered to be important for the antidiuretic activity, and also some tendency to adopt a 5,6 ß-turn. In desmopressin, in contrast to the native vasopressin, deamino-vasopressin and [D-Arg8]-vasopressin (DAVP), the Arg8 side chain, crucial for the pressor and antidiuretic activities, is very well exposed for interaction with the receptor, whereas Gly9, crucial for the pressor and uterotonic activities, is situated together with the C-terminal amide group very close to the tocin ring. The arrangements of the Gln4 and Asn5 side chains, being crucial for OT activity, also differ in desmopressin as compared to those of AVP, dAVP and DAVP. These differences in arrangement of the important for activities side chains are likely to explain extremely potent and selective antidiuretic activities of desmopressin. © 2016 Wiley Periodicals, Inc. Biopolymers (Pept Sci) 106: 245-259, 2016.

Intra- and intersubunit changes accompanying thermal activation of the HtrA2(Omi) protease homotrimer.

HtrA2(Omi) protease is involved in the maintenance of mitochondrial homeostasis and stimulation of apoptosis as well as in development of cancer and neurodegenerative disorders. The protein is a homotrimer whose subunits comprise serine protease domain (PD) and PDZ regulatory domain. In the basal, inactive state, a tight interdomain interface limits access both to the PDZ peptide (carboxylate) binding site and to the PD catalytic center. The molecular mechanism of activation is not well understood. To further the knowledge of HtrA2 thermal activation we monitored the dynamics of the PDZ-PD interactions during temperature increase using tryptophan-induced quenching (TrIQ) method. The TrIQ results suggested that during activation the PDZ domain changed its position versus PD inside a subunit, including a prominent change affecting the L3 regulatory loop of PD, and also changed its interactions with the PD of the adjacent subunit (PD*), specifically with its L1* regulatory loop containing the active site serine. The α5 helix of PDZ was involved in both, the intra- and intersubunit changes of interactions and thus seems to play an important role in HtrA2 activation. The amino acid substitutions designed to decrease the PDZ interactions with the PD or PD* promoted protease activity at a wide range of temperatures, which supports the conclusions based on the TrIQ analysis. The model presented in this work describes PDZ movement in relation to PD and PD*, resulting in an increased access to the peptide binding and active sites, and conformational changes of the L3 and L1* loops.

Arginine-, D-arginine-vasopressin, and their inverso analogues in micellar and liposomic models of cell membrane: CD, NMR, and molecular dynamics studies.

We describe the synthesis, pharmacological properties, and structures of antidiuretic agonists, arginine vasopressin (AVP) and [D-Arg(8)]-vasopressin (DAVP), and their inverso analogues. The structures of the peptides are studied based on micellar and liposomic models of cell membranes using CD spectroscopy. Additionally, three-dimensional structures in mixed anionic-zwitterionic micelles are obtained using NMR spectroscopy and molecular dynamics simulations. NMR data have shown that AVP and DAVP tend to adopt typical of vasopressin-like peptides ß-turns: in the 2-5 and 3-6 fragments. The inverso-analogues also adopt ß-turns in the 3-6 fragments. For this reason, their inactivity seems to be due to the difference in side chains orientations of Tyr(2), Phe(3), and Arg(8), important for interactions with the receptors. Again, the potent antidiuretic activity of DAVP can be explained by CD data suggesting differences in mutual arrangement of the aromatic side chains of Tyr(2) and Phe(3) in this peptide in liposomes rather than of native AVP. In the presence of liposomes, the smallest conformational changes of the peptides are noticed with DPPC and the largest with DPPG liposomes. This suggests that electrostatic interactions are crucial for the peptide-membrane interactions. We obtained similar, probably active, conformations of the antidiuretic agonists in the mixed DPC/SDS micelles (5:1) and in the mixed DPPC/DPPG (7:3) liposomes. Thus it can be speculated that the anionic-zwitterionic liposomes as well as the anionic-zwitterionic micelles, mimicking the eukaryotic cell membrane environment, partially restrict conformational freedom of the peptides and probably induce conformations resembling those of biologically relevant ones.

Structural studies of the C-terminal 19-peptide of serum amyloid A and its Pro → Ala variants interacting with human cystatin C.

Serum amyloid A (SAA) is a multifunctional acute-phase protein whose concentration in serum increases markedly following a number of chronic inflammatory and neoplastic diseases. Prolonged high SAA level may give rise to reactive systemic amyloid A (AA) amyloidosis, where the N-terminal segment of SAA is deposited as amyloid fibrils. Besides, recently, well-documented association of SAA with high-density lipoprotein or glycosaminoglycans, in particular heparin/heparin sulfate (HS), and specific interaction between SAA and human cystatin C (hCC), the ubiquitous inhibitor of cysteine proteases, was proved. Using a combination of selective proteolytic excision and high-resolution mass spectrometry, a hCC binding site in the SAA sequence was determined as SAA(86-104). The role of this SAA C-terminal fragment as a ligand-binding locus is still not clear. It was postulated important in native SAA folding and in pathogenesis of AA amyloidosis. In the search of conformational details of this SAA fragment, we did its structure and affinity studies, including its selected double/triple Pro → Ala variants. Our results clearly show that the SAA(86-104) 19-peptide has rather unordered structure with bends in its C-terminal part, which is consistent with the previous results relating to the whole protein. The results of affinity chromatography, fluorescent ELISA-like test, CD and NMR studies point to an importance of proline residues on structure of SAA(86-104). Conformational details of SAA fragment, responsible for hCC binding, may help to understand the objective of hCC-SAA complex formation and its importance for pathogenesis of reactive amyloid A amyloidosis.

The LA loop as an important regulatory element of the HtrA (DegP) protease from Escherichia coli: structural and functional studies.

Bacterial HtrAs are serine proteases engaged in extracytoplasmic protein quality control and are required for the virulence of several pathogenic species. The proteolytic activity of HtrA (DegP) from Escherichia coli, a model prokaryotic HtrA, is stimulated by stressful conditions; the regulation of this process is mediated by the LA, LD, L1, L2, and L3 loops. The precise mechanism of action of the LA loop is not known due to a lack of data concerning its three-dimensional structure as well as its mode of interaction with other regulatory elements. To address these issues we generated a theoretical model of the three-dimensional structure of the LA loop as per the resting state of HtrA and subsequently verified its correctness experimentally. We identified intra- and intersubunit contacts that formed with the LA loops; these played an important role in maintaining HtrA in its inactive conformation. The most significant proved to be the hydrophobic interactions connecting the LA loops of the hexamer and polar contacts between the LA' (the LA loop on an opposite subunit) and L1 loops on opposite subunits. Disturbance of these interactions caused the stimulation of HtrA proteolytic activity. We also demonstrated that LA loops contribute to the preservation of the integrity of the HtrA oligomer and to the stability of the monomer. The model presented in this work explains the regulatory role of the LA loop well; it should also be applicable to numerous Enterobacteriaceae pathogenic species as the amino acid sequences of the members of this bacterial family are highly conserved.

Conformational studies of neurohypophyseal hormones analogues with glycoconjugates by NMR spectroscopy.

Two glycosylated peptides have been studied using NMR spectroscopy supported by molecular modeling. Peptide I is an oxytocin (OT) analogue in which glutamine 4 was replaced by serine with attached α-d-mannose through the oxygen ß atom, whereas peptide II is a lysine-vasopressin (LVP) analogue with lysine 8 side chain modified by the attachment of glucuronic acid through an amide bond. Both peptides exhibit very weak uterotonic effect and are less susceptible to proteolytic degradation than the mother hormones. Additionally, peptide II reveals very weak pressor and antidiuretic activities. Our results have shown that the conformational preferences of glycosylated analogues are highly similar to those of their respective mother hormones. OT glycosylated analogue (I) exhibits a 3,4 ß-turn characteristic of OT-like peptides, and vasopressin-glycosylated analogue (II) exhibits ß-turns typical of vasopressin-like peptides. Therefore, the lack of binding of the glycosylated analogues to the receptors can be attributed to a steric interference between the carbohydrate moieties and the receptors. We also consider this to be the reason of the very low activity of the analyzed glycopeptides. We expect that results from these studies will be helpful in designing new OT-like and vasopressin-like drugs.

Temperature-induced changes of HtrA2(Omi) protease activity and structure.

HtrA2(Omi), belonging to the high-temperature requirement A (HtrA) family of stress proteins, is involved in the maintenance of mitochondrial homeostasis and in the stimulation of apoptosis, as well as in cancer and neurodegenerative disorders. The protein comprises a serine protease domain and a postsynaptic density of 95 kDa, disk large, and zonula occludens 1 (PDZ) regulatory domain and functions both as a protease and a chaperone. Based on the crystal structure of the HtrA2 inactive trimer, it has been proposed that PDZ domains restrict substrate access to the protease domain and that during protease activation there is a significant conformational change at the PDZ-protease interface, which removes the inhibitory effect of PDZ from the active site. The crystal structure of the HtrA2 active form is not available yet. HtrA2 activity markedly increases with temperature. To understand the molecular basis of this increase in activity, we monitored the temperature-induced structural changes using a set of single-Trp HtrA2 mutants with Trps located at the PDZ-protease interface. The accessibility of each Trp to aqueous medium was assessed by fluorescence quenching, and these results, in combination with mean fluorescence lifetimes and wavelength emission maxima, indicate that upon an increase in temperature the HtrA2 structure relaxes, the PDZ-protease interface becomes more exposed to the solvent, and significant conformational changes involving both domains occur at and above 30 °C. This conclusion correlates well with temperature-dependent changes of HtrA2 proteolytic activity and the effect of amino acid substitutions (V226K and R432L) located at the domain interface, on HtrA2 activity. Our results experimentally support the model of HtrA2 activation and provide an insight into the mechanism of temperature-induced changes in HtrA2 structure.

Highly potent antidiuretic antagonists: conformational studies of vasopressin analogues modified with 1-naphthylalanine enantiomers at position 2.

In this paper, we investigated the structure-activity relationship of two vasopressin analogues, [Cpa(1),(L-1-Nal)(2) ]AVP (I) and [Cpa(1),(D-Nal)(2) ]AVP (II) by NMR spectroscopy and molecular modeling. Both peptides exhibit antioxytocic and antipressor potency. Inversion of configuration of the residue at position 2 converted a weak antidiuretic agonist (peptide I) into a highly potent antidiuretic antagonist (peptide II). For this reason, the purpose of our study was to explain the causes of different interactions of the analogues with V(2) receptors. The results have shown that both analogues display the tendency to adopt ß-turns in the 1-4 and 2-5 fragments, which is characteristic of OT and V(1a) receptors antagonists. In addition, the [Cpa(1),(L-1-Nal)(2) ]AVP (I) shows the propensity to assume ß-turn at position 7,8, which is believed to enhance antidiuretic activity, although not being crucial for its appearance. Moreover, the C-terminal amide group seems to be crucial for signal transduction. Its high accessibility in [Cpa(1),(L-1-Nal)(2) ]AVP (I) in contrast to [Cpa(1),(D-1-Nal)(2) ]AVP (II), probably results in V(2) receptor activation.

Ensemble fits of restrained peptides' conformational equilibria to NMR data. Dependence on force fields: AMBER/8 ff03 versus ECEPP/3.

Two variants of NMR-based conformational analyses of flexible peptides are compared using two examples meeting the formula Tyr-D-Daa-Phe-Daa-NH2 (Daa=diamino acid): 1 combining D-Dab² (α,γ-diaminobutyryl) with Lys4, and 2 -D-Dap² (α,ß-diaminopropionyl) with Orn4. The ω-amino groups of D-Daa² and Daa4 are coupled with C=O into the urea, restraining 1 and 2 with 16- and 14-membered rings and leading to potent and impotent µ/δ opioid peptides, respectively. To the current task, we took from an earlier work (Filip et al, J. Pept. Sci. 11 (2005) 347-352) the NMR NOE- and J-data in H2O/D2O; and the selection of the ensembles of 1 and 2, 822 and 788 conformational families, respectively, obtained by using the EDMC/ECEPP3 method. Here, we generated ensembles of 1 and 2 using AMBER molecular dynamics in explicit water to eventually selected 686 and 761 conformers for 1 and 2, respectively. We did numbers of fits for both types of the conformational ensembles of 1 and 2 to their NOE- and J-data using a common method i.e. maximum entropy approach (Groth et al, J. Biomol. NMR 15 (1999) 315-330). Both types of the well structurally diversified ensembles fit to quite different equilibria in regressions to common experimental NOE- and J-restraints using maximum entropy principle, which is a disappointing message. Intriguing is startlingly small standard deviation in J-couplings: σ(JNHαH) ≈ 0.01 Hz for LES-MD/AMBER ensemble, contrary to σ(JNHαH) = 0.8 - 1.1 Hz for the EDMC/ECEPP ensemble, over the wide range of entropy, i.e. relatively insensitive to it. A similar feature is not the case when comparing σ(NOE) in both methods. Hence, at minute entropy contributions, it follows that J does or does not transpose "overfitted" into the final σ(J) in the AMBER or ECEPP ensemble, respectively. Could this be an effect of softness of the AMBER flexible-valence force field compared to ECEPP rigid-geometry, and its effect on ensemble sampling? We do not know an answer.

Synthesis, biological activity and solution structure of new analogues of the antimicrobial Gramicidin S.

Gramicidin S (GS) is a cyclo-decapeptide antibiotic isolated from Bacillus brevis. The structural studies have shown that GS forms a two-stranded antiparallel ß-sheet imposed by two II' ß-turns. Despite its wide Gram+ and Gram- antimicrobial spectrum, GS is useless in therapy because of its high hemotoxicity in humans. It was found, however, that the analogues of GS-14 (GS with 14 amino acid residues) attained a better antimicrobial selectivity when their amphipatic moments were perturbed. In this study, we report effects of similar perturbations imposed on GS cyclo-decapeptide analogues. Having solved their structures by NMR/molecular dynamics and having tested their activities/selectivities, we have concluded that the idea of perturbation of the amphipatic moment does not work for GS-10_0 analogues. An innovative approach to the synthesis of head-to-tail cyclopeptides was used.

NMR studies of new arginine vasopressin analogs modified with alpha-2-indanylglycine enantiomers at position 2 bound to sodium dodecyl sulfate micelles.

In this paper, we use NMR spectroscopy and molecular modeling to examine four new vasopressin analogs modified with alpha-2-indanylglycine (Igl) at position 2, [L-Igl(2)]AVP (I), [D-Igl(2)]AVP (II), [Mpa(1),L-Igl(2)]AVP (III) and [Mpa(1),D-Igl(2)]AVP (IV), embedded in a sodium dodecyl sulfate (SDS) micelle. All the analogs display antiuterotonic activity. In addition, the analogs with D-Igl reveal antipressor properties. Each analog exhibits the tendency to adopt beta-turns at positions 2, 3 and/or 3, 4, which is characteristic of oxytocin-like peptides. Mutual arrangement of aromatic residues at positions 2 and 3 has been found to be crucial for binding antagonists with the OT and V(1a) receptors. The orientation of the Gln(4) side chain seems to be important for the V(1a) receptor affinity. In each of the peptides studied, the Gln(4) side chain is folded back over the ring moiety. However, it lies on the opposite face of the tocin moiety in analogs with L and D enantiomers of Igl.

The role of the L2 loop in the regulation and maintaining the proteolytic activity of HtrA (DegP) protein from Escherichia coli.

The aim of this study was to characterize the role of particular elements of the regulatory loop L2 in the activation process and maintaining the proteolytic activity of HtrA (DegP) from Escherichia coli. We measured the effects of various mutations introduced to the L2 loop's region (residues 228-238) on the stability of HtrA molecule and its proteolytic activity. We demonstrated that most mutations affected the activity of HtrA. In the case of the following substitutions: L229N, N235I, I238N, the proteolytic activity was undetectable. Thus, the majority of interactions mediated by the studied amino-acid residues seem to play important role in maintaining the active conformation. Formation of contacts between the apical parts (residues 231-234) of the L2 loops within the HtrA trimer, in particular the residues D232, was shown to play a crucial role in the activation process of HtrA. Stabilization of these intermolecular interactions by substitution of D232 with valine caused a stimulation of proteolytic activity whereas deletion of this region abolished the activity. Since the pathogenic E. coli strains require active HtrA for virulence, the apical part of L2 is of particular interest in terms of structure-based drug design for treatment E. coli infections.

Temperature-induced conformational changes within the regulatory loops L1-L2-LA of the HtrA heat-shock protease from Escherichia coli.

The present investigation was undertaken to characterize mechanism of thermal activation of serine protease HtrA (DegP) from Escherichia coli. We monitored the temperature-induced structural changes within the regulatory loops L1, L2 and LA using a set of single-Trp HtrA mutants. The accessibility of each Trp residue to aqueous medium at temperature range 25-45 degrees C was assessed by steady-state fluorescence quenching using acrylamide and these results in combination with mean fluorescence lifetimes (tau) and wavelength emission maxima (lambda(em)max) were correlated with the induction of the HtrA proteolytic activity. Generally the temperature shift caused better exposure of Trps to the quencher; although, each of the loops was affected differently. The LA loop seemed to be the most prone to temperature-induced conformational changes and a significant opening of its structure was observed even at the lowest temperatures tested (25-30 degrees C). To the contrary, the L1 loop, containing the active site serine, remained relatively unchanged up to 40 degrees C. The L2 loop was the most exposed element and showed the most pronounced changes at temperatures exceeding 35 degrees C. Summing up, the HtrA structure appears to open gradually, parallel to the gradual increase of its proteolytic activity.

Molecular dynamics study of the internal water molecules in vasopressin and oxytocin receptors.

The role of the internal water molecules in vasopressin and oxytocin receptors has been investigated via molecular dynamics simulations in hydrated membrane model. Several water molecules have been identified within the binding pockets of receptors, where they interact with the conserved residues. In all unliganded receptors, the water molecules bound to the highly conserved D2.50 cluster have been observed. It has been proposed which water molecules may significantly contribute to the stability of overall receptor structure. In receptor-ligand complexes the water molecules are involved in the receptor-ligand interactions by forming water-mediated hydrogen bonds at their contact surface.

Conformation-activity relationships of cyclo-constrained micro/delta opioid agonists derived from the N-terminal tetrapeptide segment of dermorphin/deltorphin.

The N-terminal tetrapeptide segments of dermorphin (Tyr-D-Ala-Phe-Gly-Tyr-Pro-Ser-NH(2)) and deltorphin (Tyr-D-Ala-Phe-Asp/Glu-Val-Val-Gly-NH(2)) are agonists at the opioid receptors micro and delta, respectively. [D-Arg(2), Lys(4)]-dermorphin-(1-4) amide (Tyr-D-Arg-Phe-Lys-NH(2), DALDA) and [Dmt(1)]DALDA (where Dmt is 2',6'-dimethyltyrosine) are among the most potent and selective micro-agonists reported to date, both in vitro (having picomolar micro receptor affinity) and in vivo. In this communication, conformation-activity studies of the following four cyclic analogs of DALDA are presented and discussed: the lead peptide S(2),S(4)-cyclo (Tyr-D-Cys-Phe-Cys-NH(2)), constrained by means of an S(4.2)--S(4.4) disulfide between Cys(2) and Cys(4); its two cis and trans C(4.2)--C(4.4)-olefinic dicarba analogs, and the product of saturation of them both. They are potent nonselective or moderately micro-selective opioid agonists in vitro.They have been synthesized and tested earlier [Berezowska I, Chung NN, Lemieux C, Wilkes BC, and Schiller PW, Acta Biochim Polon 53, 2006, 73-76]. We have studied their conformations using NMR and molecular dynamics. With major conformational constraints imposed by the 11-membered ring spanning residues 2-4, they show well defined conformations of this ring, while the exocylic Tyr(1) and Phe(3) side chains still have significant conformational freedom. The more active and selective micro versus delta disulfide and saturated dicarba agonists seem to have in common: (i) their ring structures more flexible than those of the other two and (ii) their ring structures similar to each other and more diverse than those in the other two. Given this and the small size of the peptides having confirmed bioactivity profiles, there is a chance that their conformations determined in solution approach receptor-bound conformations.

Deltorphin analogs restricted via a urea bridge: structure and opioid activity.

Eight cyclic heptapeptides related to the full sequence of deltorphin have been synthesized. The synthesis of linear peptides containing diamino acid residues in positions 2 and 4 was carried out on a 4-methylbenzhydrylamine resin. Depending on protection procedures, the N-protected peptide-resins or N-protected peptide amides with free amino groups in the side chains were obtained, which were subsequently treated with bis-(4-nitrophenyl)carbonate to form a urea unit. Opioid activities of the peptides were determined in the guinea pig ileum (GPI) and mouse vas deferens (MVD) assays. Several compounds showed high delta opioid agonist potency and high selectivity for delta receptors. The results were compared with those obtained earlier for respective 1-4 deltorphin analogs. The conformations of these peptides have been studied using 2D-NMR in H2O/D2O and molecular dynamics. We observed that the backbone rings had well defined conformations, while the Tyr and Phe side chains and the C-terminal tail had significant conformational freedom. The bioassay data and conformational parameters of these peptides were compared with those of previously described, corresponding 1-4 deltorphin analogs. This comparison permitted an assessment of the role of the C-terminal peptide segment in defining the conformation and receptor interaction of the N-terminal portion and provided insight into the relationship between the putative bioactive conformations and bioactivity.

An enormously active and selective azapeptide inhibitor of cathepsin B.

The peptidomimetic Z-Arg-Leu-Arg-Agly-Ile-Val-OMe (where Agly means alpha-aza-glycyl, -NHNHCO-) is the strongest (K(i) = 480 pM) and the most selective inhibitor of cathepsin B to date, being approximately 2310 times as active to cathepsin B as to cathepsin K. In this paper we introduce the peptide and seek to rationalize its structure-activity relationships using molecular dynamics (MD) and NMR. It is shown that the -Agly-moiety restrains the peptide backbone to a bent shape, contrary to its parent peptide (with Gly in position 4), having its backbone extended and flexible. This fold is maintained in the plug covalently bound to the cathepsin B Cys29, in compliance with similar bends already observed in two other azapeptides attached to the active sites of cathepsin B. The MD simulation of the Z-Arg-Leu-Arg-Agly approximately cathepsin B complex suggests that, contrary to other potent inhibitors of cathepsin B, the current double Arg(1)/Arg(3) inhibitor, while maintaining the fold is able to form a unique ion cluster involving both Arg residues on the inhibitor part and two acidic Glu171 and Glu245 on the cathepsin B part, thus enhancing the affinity and subsequently the inhibiting power and selectivity of Z-Arg-Leu-Arg-Agly-Ile-Val-OMe to the observed extreme extent.